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ABSTRACT 

Swift J1644+57 is an unusual transient event, likely powered by the tidal disruption 
of a star by a massive black hole. There are multiple short timescales X-ray flares 
were seen over a span of several days. We propose that these flares could be produced 
by internal shocks. In the internal shock model, the forward and reverse shocks are 
produced by collisions between relativistic shells ejected from central engine. The 
synchrotron emission from the forward and reverse shocks could dominate at two quite 
different energy bands under some conditions, the relativistic reverse shock dominates 
the X-ray emission and the Newtonian forward shock dominates the infrared and 
optical emission. We show that the spectral energy distribution of Swift J1644-f57 
could be explained by internal shock model. 
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1 INTRODUCTION 

Swift J1644+57 was triggered by the Swift/BAT on 28 
March 2011 (Cummings et al. 2011). Swift J1644-I-57 was 
initially discovered as a long-duration gamma-ray burst 
(GRB 110328A) by the Swift satellite, but the light curve 
soon showed that it was quite different. It remained bright 
and highly variable for a long period, and re-trigger the BAT 
three times over the next 48 hours (Sakamoto et al. 2011). 
The isotropic X-ray luminosity of Swift J1644-I-57 ranges 
from 10*^ — 4 X 10** erg s~^, and the total isotropic en- 
ergy is about 3 x 10^^ erg during the first 30 days after the 
BAT trigger (Burrows et al. 2011). From the strong emis- 
sion lines of hydrogen and oxygen, Levan et al. (2011) esti- 
mate the redshift of Swift J1644+57 is z = 0.35. From the 
astrometric observation of the X-ray, optical, infrared, and 
radio transient with the light-centroid of the host galaxy, it 
is found that the position of this source is consistent with 
arising in the nucleus of the host galaxy (Bloom et al. 2011; 
Zauderer et al. 2011). 

The X-ray hght curve of Swift J1644-I-57 exhibits re- 
peated extremely short timescale flares. The flares have rise- 
times as short as 100 s (Burrows et al. 2011). These flares 
are similar as the flares discovered in the GRB afterglow 
(Burrows et al. 2005), which may indicate the same origin 
of them. The internal shock model can produce the X-ray 
flares observed in GRB afterglows (Burrows et al. 2005; Fan 
& Wei 2005; Zhang et al. 2006; Yu & Dai 2009). 

After the Swift J1644-I-57 was discovered, several mod- 
els were proposed to explain it, most concentrating on a 



picture that a main sequence star was tidally disrupted by 
passing too close to a 10® — 10^ Mq black hole (Bloom et 
al. 2011; Burrows et al. 2011; Cannizzo et al. 2011; Socrates 
2011; Shao et al. 2011). Krohk & Phan (2011) suggest that 
this event may be produced by a white dwarf tidally dis- 
rupted by a 10* M0 black hole. The process is as follow: a 
star is disrupted as it passes near a supermassive black hole, 
and much of its mass is distributed into an accretion disk 
around the black hole. A powerful jet is then launched. In 
these models, the X-ray emission is thought to be produced 
by external inverse Compton (EIC) (Bloom et al. 2011) or 
synchrotron emission (Burrows et al. 2011). But on the high 
frequency side, the Fermi LAT (Campana et al. 2011) and 
VERITAS upper limits (Aliu et al. 2011) require that the 
synchrotron self-Compton (SSC) component is suppressed 
by 7 — 7 pair production. The soft photons of 7 — 7 pair 
production are thought to be generated from the thermal 
emission of the accretion disk or the disk outflow. In the 
SSC model, soft photons originated from the thermal emis- 
sion of the accretion disk may not provide an efficient source 
for the 7 — 7 production. Because the condition of 7 — 7 pro- 
duction is ExE^i{l — cos9) 2(meC^)^, where 9 is the angle 
between the directions of soft seed photon and high-energy 
photon. Only a fraction of high-energy emission can be ab- 
sorbed by soft photons. So the soft photons from the disk 
outflow may be a better candidate (Strubbe & Quataert 
2009). In the synchrotron emission model, the jet must have 
a strong magnetic held (Foynting-flux-dominated) and has 
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ongoing in situ acceleration of electrons (Aliu et al. 2011; 
Burrows et al. 2011). 

In this letter, we use the internal shock model to explain 
the X-ray flares of Swift J1644-I-57. The internal shock pro- 
duces the prompt emission of GRB in the standard fireball 
model (Rees & Meszaros 1994; Paczyriski & Xu 1994). The 
internal shock model also the leading model of X-ray flares 
in GRBs, the external shock model is very hard to account 
for the X-ray flares (Burrows et al. 2005; Fan & Wei 2005). 
The central engine of this event may be formed as follow. 
When a supermassive black hole tidal disrupts a star, a disk 
is formed. The magnetic field could be produced by disk 
instability. The disk can then anchor and amplify the seed 
magnetic field to a strong ordered poloidal field , which in 
turn threads the black hole with debris material in the inner 
region of the disk. A large amount of the rotational energy 
of the black hole can be extracted via the Blandford-Znajek 
(BZ) process, which creates a jet along the rotation axis 
(Blandford & Znajek 1977). The magnetic field lines will 
break the disk into blobs, so many shells could be ejected 
(Cheng & Lu 2001). When fast shell catches up with early 
slow shell, internal shock is generated. Other models of cen- 
tral engine also discussed, such as the episodic accretion onto 
a central object due to a chopped accretion disk (Perna et 
al. 2006) , or episodic accretion due to a modulation of the 
accretion fiow by a magnetic barrier (Proga & Zhang 2006). 

The structure of this letter is as follow. In next section, 
we describe the dynamics of internal shock arising from a 
collision between two shells and the synchrotron radiation 
of the shocked electrons. In Section 3, we apply the model to 
the Swift J1644+57. Finally, a summary is given in Section 
4. 



2 THE INTERNAL SHOCK MODEL 

The internal shock model has been extensively discussed in 
literature (Rees & Meszaros 1994; Paczynski & Xu 1994; Yu 
& Dai 2009; Yu, Wang & Dai 2009). We give a brief descrip- 
tion of our model as follow. Shells with different Lorentz 
factors and densities are ejected by the central black hole. 
Collisions of a pair of ejecta can produce different intensities 
of X-rays. For example two shells with similar Lorentz fac- 
tor and density will produce a weak flare whereas two shells 
with large differences can produce a strong flare. Since the 
collision frequency of these pairs (internal shocks) should 
be very high, it should result in rapid variable intensities. 
Simultaneously some earlier ejected fast moving shells can 
already reach the interstellar medium(ISM) and produce the 
external shock there (Sari & Piran 1995). Since the earlier 
ejected shells are less and hence the radiation results from 
external shock should be weak in the beginning. However 
after collisions of pairs they can merge and move toward 
the ISM and provide more energy into the external shock. 
Therefore the radiation intensity due to the external shock 
should gradually increase. We should note that since the in- 
jected energy provided in this way is in a discrete manner, 
therefore the flux will increase substantially but gradually 
decrease back to the original light curve. This phenomenon is 
similar to that observed in GRBs known as "re-brightening" 
effect (Zhang & Meszaros 2002; Huang et al. 2006). In Fig- 
ure [1] we provide a schematic illustration of our model. 




Figure 1. Schematic illustration of the model. The dotted line 
represents the intensity produced by internal shocks. The dashed 
line represents the intensity produced by the external shock. In 
this figure we argue that large short-term fluctuations can still 
occur due to collisions of few later ejecta when the jet still exists. 



2.1 Shock dynamics 

At a time of t,4, the central engine ejects a shell denoted by 
A with bulk Lorentz factor 7a with isotropic kinetic energy 
luminosity La - Some time {5t) later, another shell B with 73 
and Lb is assumed to be ejected. In order to let shell B catch 
up and collide with the shell A, 7b > 7a is required. At the 
radius J?coi = l3Al3Bc5t/ip{z){PB — /3a), a collision between 
A and B takes place. For (7A, 7b) 3> 1, The collision radius 
is (Yu & Dai 2009) 



2'ylc5t 



(1 - (7a/7b)2)V'(^) 



(1) 



where tp{z) — 1 + z. 

After the collision, a forward shock and a reverse shock 
are produced. The system is separated into four regions by 
the two shocks and a contact discontinuity surface: (1) un- 
shocked shell A, (2) shocked shell A, (3) shocked shell B, 
and (4) unshocked shell B, bulk Lorentz factors of which 
are 71 = 7a, 72 = 73 = 7, and 74 = 7b. Two relative 
Lorentz factors of the shocked regions relative to unshocked 
regions 1 and 4, can be calculated by 
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According to Blandford & McKee (1976), the internal energy 
densities of the two shocked regions are 62 — (721 — l)(472i -f- 
3)nimpC^ and 63 = (734 — 1)(4734 -I- 3)n4mpC^, where ni = 
La /47r_R^oi 7a "T'pC^ and n4 = Ls/inRl^i^'^mpC^ . The me- 
chanical equilibrium between the two shocked regions re- 
quires 62 = 63, so 
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where Li = La and L4 = Lb. We can calculate the values of 
7, 721 and 734 from equations ([2]) and Q after the param- 
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eters of shells are given. In four limit cases, these equations 
can be solved analytically (Yu & Dai 2009). For 74 ^ 71, 
(1) if Li/Li > (1/7) (74/71)", we have 721 = 74/271 > 1, 
734 — 1 ^ I'l/Tfll and 7 — 74(1 — -\/2f)i which means the 
forward shock is relativistic and the reverse shock is New- 
tonian; (2) if 16 ^ L4/L1 <C (1/16) (74/71)", we can obtain 

_fl/4 1/2 /„ 1/2 , 1/2 /r, rl/4 1/2 , J 

721 =/ ^ 74 /27i > 1, 734 = 74 /2/ ^ 7i > 1 and 
7 = /"'^''"7i'^^7y^i so both the two shocks are relativistic; (3) 
if I/4/ii < 7, we get 721 - 1 ^ f^lHll = i, 734 = 74/271 
and 7 = 71(1 + v^), so the forward shock is Newtonian 
and the reverse shock is relativistic. Finally, (4) for 74 « 71, 
both the two shocks are Newtonian. Since 71, 74, and / are 
unchanged with the moving of the shells, the values of 7 ,721 
and 734 are constant before the shocks cross the shells (Yu 
& Dai 2009). 
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2.2 Synchrotron emission from forward and 
reverse shocks 

Following Dai & Lu (2002), the total number of the electrons 
swept-up by the forward and reverse shocks during a period 
of St can be expressed by Ne,2 = 2y/2S^LA5t/ {ili(z)^\mpC?^ 
and A^e,3 = LsSt/ {^{z)'y4,mp(?^ , respectively (Yu, Wang & 
Dai 2009). 

The forward and reverse shocks can accelerate par- 
ticles to high energies. Following Sari et al. (1998), we 
assume that the energies of the hot electrons and mag- 
netic fields are fractions and es of the total internal 
energy, respectively. Thus, the strength of the magnetic 
fields is Bi = (SvreB^iei)^''^ , i = 2,3. We assume a 
power-law distribution of the shock-accelerated electrons, 
due.jd'^e. oc 7^^ for 7e ^ 7e,m (Sari et al. 1998). The random 
Lorentz factor of electrons in regions 2 or 3 is determined 



by 7e 



(p-2) 



(r — 1), where F equals to 721 or 



734. In both shocked regions, the hot electrons with energies 
above 7e,c,i?TieC^ lose most of their energies during a cool- 
ing time 5t, where the cooling Lorentz factor is determined 
by 7e,c,i = 67rmeCi/'(z)/ (^arBf'yStj . The two characteristic 
frequencies and a peak flux density are (Sari et al . 1998; 
Wijers & Galama 1999) 
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= is the luminosity 

/n„{i+z')^+nA 

distance of the source and ^{p) is a function of p, for p = 2.2, 
$(p) « 0.6 (Wijers & Galama 1999). In the calculation, we 
use Q.M = 0.3, Q.A = 0.7 and i/o=70 km s"^ Mpc"\ is the 
electron charge and ctt is the Thomson cross section. The 
synchrotron spectrum can be written as (Sari et al. 1998) 
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Figure 2. The broadband spectral energy distribution of Swift 
J1644-I-57 at 2.9 day after the BAT trigger. The points show the 
observation data, which are taken from Bloom et al. (2011). The 
two solid lines represent the unabsorbed spectrum of the reverse 
and forward shocks, which are generated by the internal shock. 
The dashed line shows the spectrum of external shock. In order to 
fit the spectrum, a moderate extinction {Ay = 3 — 5) is required. 



v <i'i; 

(5) 



where vi = vamivrn,i,Vc.,i), Vh ~ max{um.,i,Vc,i 
for Uc,i < Vm.i and g = p for Uc,i > v-n 



I, and g = 2 



3 IMPLICATION FOR THE SWIFT J1644+57 

There are two peaks in the spectrum of Swift J1644+57, 
far-infrared (FIR) and hard X-ray peaks. In order to fit the 
spectrum, we focus on the case (3) of internal shock model in 
section 2.1, in which the reverse shock is relativistic and the 
forward shock is Newtonian. In the rest of the paper we de- 
note Q — 10^ Qx- For illustration purpose, we set L4 = Li = 
L = lO^^ Oerg s-\ 74 = 1000, 71 = 10, £^,2 = ee,3 = = 0.5 
and e_B.2 ~ es,3 = e_B = 0.1. As shown in Cheng & Lu 
(2001), the Lorentz factor of shell can be up to 1000. 
It is reasonable to adopt 74 = 1000. According to obser- 
vation, we use 5t ~ 100 s, the variability timescale of flare 
(Bloom et al. 2011; Burrows et al. 2011). The collision ra- 
dius is _Rcoi ~ 2'y1c5t/ip{z) ~ 5 x lO^'^cm, which is consistent 
with the X-ray emission radius determined from observation 
(Bloom et al. 2011). The Lorentz factor of merged shell is 
7 ~ 14. Using equation @, we can obtain the following 
expressions for the reverse shock 



Vc,3 
,max,3 



1.2 X 10^* Hz 



2 rl/2 1/2 r,-l -4 
-0.374,3-t-47 ^B-l°^2 7l,l> 



2.2 X 10'^ Hz L-^^^^eg%l5t2'y!,i, 

no T r3/2 1/2 -2 -1 ■-2 

0.9 mjy L^; £3 _i7i,i74,3az„27.7- 
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For the forward shock, we obtain 

o c w loll TU- 2 rl/2 1/2 j,-l -2 

iym,2 ~ 3.5 X 10 tiz e^^_o,3L^'^ e^ _-^dt2 7i_i, 



l^c,2 

,max,2 



2.2 X 10^^ Hz L'l^''^eg^/_X5t2-il 

, r T r3/2 1/2 -3 ,-2 

15 mJy L^y <;fl,_i7i,i"i,27.7- 



(7) 
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Therefore, the resulting synchrotron photons emitted by the 
two shocks are expected to peak at two different energy 
bands and thus two distinct spectral component^. The peak 
of reverse shock spectrum will be at hard X-ray, but peak 
of the forward shock will be at FIR. The synchrotron self- 
absorption must be taken into account. In Vm,2 < Va,2 < i^c,2 
case, the synchrotron self-absorption frequency in region 2 
reads (Panaitescu & Kumar 2000) 



2/(p+A) 
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2p-2 
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In Va.i < Vc.'A < I'm, 3 casc, the synchrotron self- 
absorption frequency in region 3 can be calculated by 
(Panaitescu & Kumar 2000) 
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5qeNe,3 
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The maximum Lorentz factor is limited by the synchrotron 
losses and is given by (Cheng & Wei 1996) 



7M,i ^ {Sqe/Biar) 



1/2 



4 X lO'^Sr^''^ 



(10) 



Another mechanism to restrict the maximum energy of an 
electron is diffusion. It turns out that maximum Lorentz 
factor restricted by diffusion is much larger than that in 
equation pO|l . So The maximal synchrotron photon energy 
can be estimated (Fan & Piran 2008) 



hqeBi 
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MeV, 



(11) 



where h is the Planck constant, F equals to 721 or 734. 

The spectrum of internal shock model is shown in Fig- 
ure[2]using above parameters during a high state. The X-ray 
spectrum of Swift J1644+57 can be generated in our model. 
A moderate extinction {Av ~ 3 — 5) is required to explain 
the spectrum. This value of extinction is reasonable in this 
case, because of this event is arising in the nucleus of host 
galajcy. This value is also consistent with that determined 
by Bloom et al. (2011) and Burrows et al. (2011). Because 
of large value of synchrotron self-absorption frequency, the 
radio emission of internal shock is suppressed. From obser- 
vations, the radio emission is from larger radius comparing 
with X-ray emission (Bloom et al. 2011; Burrows et al. 2011) . 
The interaction between the first shell ejected by central en- 
gine and the ISM results in an external shock. The radio 
emission is from large radius and can be modeled by this 
external shock, similar as GRB afterglow. 

The total energy release during this initial period is 
about Eiso ~ 10^^ erg (Bloom et al. 2011). The ISM density 
is about n ~ lOcm"'^. Following Sari et al. (1998) and Bloom 
et al. (2011), we obtain the synchrotron frequencies and peak 



^ We can see that these characteristic frequencies, i.e. Vc and Urn, 
are very sensitive to the Lorentz factor. However, Kobayashi et 
al. (1997) have shown that the radiation loss is less than 10% 
of total energy, therefore there is virtually no evolution of these 
spectral parameters during the collision. 
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Figure 3. The spectral energy distribution of Swift J1644-I-57 at 
31 hours after BAT trigger. The points show the observation data, 
which arc taken from Burrows et al. (2011). The two lines repre- 
sent the unabsorbed spectrum of the reverse and forward shocks, 
which are generated by the internal shock. The parameters of two 
shocks are given in the text. 



flux of external shock as follow 



Va 



.2.0xl0^« Hz e-Ue'J^_,EU'nr, 
3.0 X 10" Hz el_,e'J^_,EU%!^,\ 

on^, inl3TT -3/2 r-,-1/2 -1,-1/2 
! 8.0 X 10 Hz Es, -2-^53 "-1 ' 

: 170 mJy e]^^_2E53ny^t~^^*dl%^^y. 



(12) 



The expression of F^,max is a little different from that of 
Sari et al. (1998), because the observer has already observed 
the edges of the jet, as discussed in Bloom et al. (2011). 
The spectrum of the external shock is shown as dashed line 
in figure [51 which can only produce a simple power law in 
X-ray region. The radio light curve is different from t^^'^'^ 
behavior observed the late X-ray light curve (Giannios & 
Metzger 2011; Metzger et al. 2011) because the radio light 
curve should be determined by the evolution of the external 
shock and has nothing to do with the accretion rate in the 
disk. 

So our model predicates that, during the high state 
(flaring), the emission of internal shock will dominate at X- 
ray band and a broken power-law spectrum is shown. During 
the low state (no flares) , there is no internal shocks and the 
emission is from the external shock. Because the typical fre- 
quencies of external shock is low at the first few days (see 
equation (|12|) '). the spectrum at X-ray band exhibits as a sin- 
gle power-law if no energy injection happens. In both states, 
the radio emission is from the external shock. 

In Figure|3l we fit the spectrum data from Burrows et al. 
(2011) detected at the same period. We adopt the parame- 
ters as follow: L4 = Li = L = lO^^'^erg s"\ 74 = 1000, 
71 = 10, 5t = 100 s, e, = 0.8, and es = 0.001. Since 
shocks produced by collisions are highly nonlinearly pro- 
cesses, therefore the microscopic parameters, i.e. te and e_g, 
can be different for different collisions. 

The X-ray flux from Swift J1644-)-57 is observed to 
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track the X-ray hardness (Bloom et al. 2011; Levan et al. 
2011). The X-ray flux and photon index exhibit a strong 
anti-correlation. This signature is a natural consequence of 
our model. In the earlier stage when the the internal shocks 
is dominated the X-ray flux is high but variable and harder. 
At later time when the external shock is dominated, the 
X-ray flux becomes lower but less fluctuating and softer. 

The durations of flares are very complicated, similar as 
the X-ray flares in GRBs. For individual flare, the duration 
can be roughly estimated as A/c, where A is the width of 
shell (Maxham & Zhang 2009). If the ejecta are coming from 
the disk around the black hole, it should have the size of disk 
~ 3rs ~ 6GM /c^ ~ 8 x IO^Mk cm, where M is the mass 
of black hole. From the minimum rise time, Bloom et al. 
(2011) and Burrows et al. (2011) have estimated Me, ~ 10. 
So the duration of individual flare should be of order of 
A/c ~ rd/c ~ 200 s. However, flares can superimpose on 
each other if shells collide near the same time. For example, 
the duration of the flare detected at 111045 s after BAT 
trigger is about 300 s, which is consistent with the rise time 
scale. But the duration of the flare at about 1.115x10^ s 
after BAT trigger with minimum rise time is longer than 
1000 s. We believe that this flare is a superposition of several 
flares. It is interesting to note that the flares detected in 
GRBs indicate that the shell width A broadens with ejected 
time. A natural broadening mechanism is shell spreading. 
After a shell enters the spreading regime, the width of the 
shell is proportional to the radius, so that if the collision 
radius is larger the duration of X-ray flare can last longer. In 
this event, we can also see that the width of flares broadens 
with ejected time. This is similar to some central engine 
models of GRBs, for example, in the fragmented disk model 
proposed by Perna et al. (2006), the clumps at larger radius 
have lower densities and tend to be more spread out so that 
the accretion time scale is longer. 



4 SUMMARY 

In this letter, we propose the internal shock model to ex- 
plain the X-ray flares of Swift J1644+57. In the internal 
shock model, collisions between a series of relativistic shells 
generate many pairs of forward and reverse shocks. The 
synchrotron emission produced by the forward and reverse 
shocks could dominate at two quite different energy bands 
if the Lorentz factors of these two types of shocks are signif- 
icantly different from each other. We show that the spectral 
energy distribution of Swift J1644-I-57 could be fitted in in- 
ternal shock model, in which the reverse shock is relativistic 
and the forward shock is Newtonian. A moderate extinc- 
tion {Av = 3 — 5) is required, this value is consistent with 
that used in Bloom et al. (2011) and Burrows et al. (2011). 
Burrows et al. (2011) showed that the high frequency spec- 
trum is produced by the synchrotron and SSC mechanisms, 
similar to poynting flux-dominated blazar jet model. The 
radio fluxes come from a larger region of the other jet. This 
model requires continuous in situ re-acceleration of electrons 
to maintain a low energy cut-ofi^ in the electron distribution 
(Aliu et al. 2011). Bloom et al. (2011) presented two mod- 
els for the spectrum: one is two-component blazar emission 
model, the other is forward shock emission from jet-ISM 
interaction plus EIC emission model. But on the high fre- 



quency side, the LAT and VERITAS upper limits require 
that the SSC component is suppressed by 7 — 7 pair pro- 
duction. The soft photons from the disk outflow may provide 
sources for the 7 — 7 production. 

The rapid rise and decline of the light curve may indi- 
cate the internal shock origin of these flares. The external 
shock is very hard to account for the X-ray flares (Burrows 
et al. 2005; Fan & Wei 2005; Zhang et al. 2006). During the 
high state, the emission of internal shock will dominate at 
X-ray band and a broken power-law spectrum is shown. Dur- 
ing the low state, there is no internal shock and the emission 
is from external shock. The spectrum at X-ray band will be 
shown as a single power-law if no energy injection happens. 
In both states, the radio emission is from the external shock. 
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